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Abstract 

By combining analytical and numerical methods we find that the solutions of 
the complete Horava theory with negative cosmological constant that satisfy the 
conditions of staticity, spherical symmetry and vanishing of the shift function are 
two kinds of geometry: (i) a wormhole-like solution with two sides joined by a 
throat and (ii) a single side with a naked singularity at the origin. We study the 
second-order effective action. We consider the case when the coupling constant of 
the (9 In N) 2 term, which is the unique deviation from general relativity in the ef¬ 
fective action, is small. At one side the wormhole acquires a kind of deformed AdS 
asymptotia and at the other side there is an asymptotic essential singularity. The 
deformation of AdS essentially means that the lapse function N diverges asymptot¬ 
ically a bit faster than AdS. This can also be interpreted as an anisotropic Lifshitz 
scaling that the solutions acquire asymptotically. 
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1 Introduction 


The nonrenormalizability of general relativity (GR) in the perturbative scheme as well as 
the issues of dark matter and dark energy open the doors for the possibility of studying 
completions of GR. In this sense the Horava proposal [1] of having an ultraviolet com¬ 
pletion of GR with a preferred foliation of space-time has been studied widely. Among 
the several versions that have been proposed, the complete nonprojectable version [1, 2] 
exhibits a good behavior in what concerns to the structure of the field equations and 
constraints [2, 3, 4, 5, 6]. For this version most of the analysis has been done outside the 
conformal point of the kinetic term, which holds at the value A = 1/3 of its corresponding 
coupling constant. Out from the conformal point, A 7^ 1/3, the theory propagates one 
extra mode in addition to the two tensorial modes of GR 4 . The decoupling of this extra 
mode at low energies faces with the so-called strong coupling problem [10, 11, 12, 13]. 
However, at the conformal point of the kinetic term, A = 1/3, the theory propagates 
exactly the same modes of GR as a consequence of additional second-class constraints 
that arise at this point [6] (see also [14]). This is a outstanding feature, since there is no 
need of any decoupling mechanism and there are no discontinuity issues at low energies. 
This suggests that theory can be smoothly connected to GR [7, 6]. Because of this we 
consider the A = 1/3 case of the nonprojectable theory as a interesting model among the 
various proposals known for quantum gravity. 

Mainly because of their potential astrophysical applications, in the modifications of 
GR it is mandatory to identify solutions representing isolated sources. A suitable sector 
to start with is the set of spherically symmetric and static configurations. Obviously, 
the most important information for these configurations comes from the large-distance 
effective theory. Thus, our main interest in this paper is to study static spherically 
symmetric solutions of the effective theory of the complete nonprojectable Horava theory. 
Since these configurations are static (and since we switch off the shift function), they do 
not depend on the constant A, thus the configurations are the same for the A = 1/3 and 
A 7^ 1/3 cases. 

The static spherically symmetric solutions with vanishing shift function of the effective 
theory without cosmological constant were found in Ref. [15] (see also [16]). Actually, 
the solutions of Ref. [15] were found in the context of the Einstein-aether theory [17]. 
The correspondence comes from the equivalence this theory has with the second-order 
effective action of the Horava theory [12, 18, 19]. An independent analysis of the same 
configurations but directly on the Horava theory was done in Ref. [20], where several 
features of the solutions were studied. Later on, we considered again the problem of the 
static spherically symmetric solutions of the complete nonprojectable Horava theory in 
Ref. [21] (with vanishing shift function). The main feature of the solutions (for positive 
mass) is that they have a wormhole-like geometry with a throat joining the two sides. 
They have an asymptotically flat side and an essential singularity at the infinite boundary 
of the other side. In Ref. [21] we recovered the same solutions of Ref. [15] but explicitly 

4 An exception (at the level of classical actions) is the model with only a J7-term in the potential, 
which is equivalent to GR for all values of A and for asymptotically flat geometries. This was shown in 
[7] and corroborated in [8, 9]. 
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on a coordinate system valid at the throat. In particular this allows to show directly that 
the physical radius r has a minimum at the throat. 

The results of our previous work [21] encourage us to consider the problem of finding 
the static spherically symmetric solutions of the Horava theory now with a cosmological 
constant turned on. Specifically, we shall consider the case of negative cosmological 
constant, under which we can get better control of the behavior of the held equations. 
As we have already mentioned, for the analysis we take the lowest-order effective action, 
which contains the cosmological constant and the terms of second order in derivatives. 
We comment that in Ref. [20] it was also studied the same effective action with negative 
cosmological constant, with focus on the asymptotic behavior of the solutions. In this 
paper we go beyond the results of [20] by performing analysis of the solutions over all 
the space with the aim of understanding their whole geometry. Our program is the 
implementation of a procedure parallel to the one of the A = 0 case we previously 
studied [21]. In particular we pursue a new radial coordinate that allows to identify the 
geometry of the solutions. We shall arrive at differential equations determining both the 
new coordinate and the metric components in terms of it. With analytical techniques 
we shall extract important properties from these equations, such as the presence of a 
minimum for the physical radius r, which indicates that the geometry has a throat. 
Since in this case the relevant equations are considerably more involved than in the 
A = 0 case, eventually we shall perform numerical integration on them. This will confirm 
the analytical results and complete the understanding of the configurations. 

It is worth stressing that, because the smaller group of symmetries and the different 
dynamics, in Horava theory one does not expect to find spherically symmetric configu¬ 
rations possessing the universality properties that the Schwarzschild solution has in GR. 
To start with, in (vacuum) Horava theory spherical symmetry does not imply staticity or 
asymptotic flatness. Therefore, in our study we impose spherical symmetry and staticity 
independently. We also explore the asymptotia of the solutions, both analytically and 
numerically, rather than imposing it. Moreover, once the general ansatz satisfying spher¬ 
ical symmetry and staticity has been put in spherical and adapted-time coordinates, it 
still has the degree of freedom of the shift function open. Specifically, the time-radial 
component is left free and it can not be gauge-fixed. Different forms of this component 
lead in general to inequivalent configurations in the context of the Horava theory. For 
the sake of simplicity we switch off this function. By doing so we obtain a def ini te kind 
of solutions, but the reader should take into account that there are other inequivalent 
solutions that are static and spherically symmetric. In particular, black holes have been 
found under such conditions in the complete nonprojectable Horava theory [22, 23]. 5 
Universal horizons which arise in solutions like those have been the subject of intense 
study, see for example [24] and references therein. Directly in the original Horava theory 
without the terms of Ref. [2], solutions with the radial component of the shift function 
turned on have also been found [25]. Other exact solutions of various versions of the 
Horava theory can be found for example in Refs. [26, 27, 28]. In particular, in Ref. [29] 
wormhole configurations were considered in the original Horava theory without the terms 

5 Since in Refs. [22, 23] the theory is formulated in a covariant way with an additional vector field, the 
extra degree of freedom arises as the radial component of the vector field rather than in the shift vector. 
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of Ref. [2], 

The study we present here can be also of interest for holography. Solutions of Horava 
theory has been used as holographic duals of Lifshitz-scaling field configurations; see, 
for example, Ref. [30]. Since we study the theory with negative cosmological constant, 
the spacetimes we study here could enter into the arena of holographic duals (at least 
numerically). 

As we have mentioned, another important aspect we shall investigate is the asymptotic 
behavior of the solutions. Under a hypothesis of polynomial divergence, which can be 
posed even before finding the global structure, we shall determine analytically whether 
the solutions are asymptotically anti-de Sitter and we shall contrast the analytical result 
with the numerical solution. We give in advance that the asymptotia we find (at one of 
the sides of the throat) is not exactly anti-de Sitter, but rather a deformation of it. At 
the other side there is an essential singularity similar to the one found in the inner side 
of the A = 0 case [15, 21]. 


2 Solutions of the field equations 

2.1 The field equations and the ansatz 

In Horava theory a given foliation of spacetime is considered of absolute physical meaning, 
such that the allowed symmetry transformations can not change the foliation [1]. The 
foliation can be defined by the existence of a global “time” function and for convenience 
the time coordinate t can be identified with this function. Then the spacetime is foliated 
by hypersurfaces of t = constant which are spacelike. In this context it is convenient 
to parameterize the spacetime metric in terms of the Arnowitt-Deser-Misner (ADM) 
variables, such that 


ds 2 = -(N 2 - NiN^dt 2 + 2 Nidx'dt + g ij dx i dx j . (2.1) 

We assume that the lapse function N depends both on time and space, hence we deal 
with the nonprojectable formulation of the theory [1], In particular, this assumption is 
fundamental to find the kind of configurations we are interested in. The most general 
kind of diffeomorphisms over the spacetime that do not change the given foliation are the 
foliation-preserving diffeomorphisms, under which the coordinates and ADM variables 
transform as 

8t = f(t), 5x z = ( l (t, x ), 

8N = ( k d k N + fN + fN , 

. . . (2.2) 
SNi = C k d k N t + N k di ( k + C j 9ij + fNi + fNi , 

$9ij C d k gij T 2 g k (idj)C, + f 9%j ■ 

Technically, the difference with the group of general diffeomorphisms used in GR is that 
the function / determining the transformation of the time is not allowed to depend on 
the spatial coordinates. That is, under the foliation-preserving diffeomorphisms the time 
coordinate transforms at most into itself, which keeps the foliation. 
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The complete nonprojectable theory was formulated in Refs. [1, 2], where the details 
of the formulation can be found. The main idea is to insert into the potential all the 
terms that preserve the gauge symmetry declared as the fundamental one, which is the 
foliation-preserving-diffeomorphisms symmetry, up to a given order in derivatives [2], 
For the renormalization of the theory it is important to include terms of, at least, sixth 
order in spatial derivatives [1]. But for solutions which are of interest for large-distance 
applications one can deal with lower-order effective actions, since the influence of the 
higher order terms on these configurations can be neglected. This procedure simplifies 
the computations enormously. When dealing with effective actions all the terms of the 
given order and below that preserve the symmetry must be included. Following the spirit 
of effective theories, each term that is separately covariant under the gauge symmetry is 
written with an independent coupling constant. For example, the kinetic part has two 
terms that are separately covariant, K t jK lJ and K 2 , where is the extrinsic curvature 
tensor for each spatial slide of the foliation, 

Kij — — 2 V(i-ZVj)), (2-3) 

and K is its trace, K = g lJ K, 3 . Therefore, a relative multiplicative constant between 
KijK and K 2 is needed, which is denoted by A. The kinetic term arises in the action as 

K, tj K ij - A K 2 . (2.4) 

In principle any value of A is valid, unlike what happens in GR where the symmetry 
of full spacetime diffeomorphisms demands A = 1 since K tJ K l] and K 2 are not sepa¬ 
rately covariant under this higher symmetry, they are only covariant when they are taken 
together. 

In the complete nonprojectable theory the value of A plays a crucial role in determining 
the number of propagating degrees of freedom. When A ^ 1/3 the theory propagates 
three degrees of freedom: the two tensorial modes that are also present in GR plus one 
extra scalar mode; see, for example, Refs. [2, 3, 4], But when A = 1/3 the theory acquires 
two extra second-class constraints that eliminate the extra scalar degree, thus at A = 1/3 
the theory propagates exactly the same degrees of freedom of GR [6]. Interestingly, at 
A = 1/3 the full kinetic term (2.4) acquires a conformal covariance [1], but in general the 
terms in the potential break this conformal symmetry. 

The most general effective action for large distances (second order in derivatives) with 
cosmological constant is [1, 2] 

S = I dtd 3 Xy/gN(K i:j K i2 - A K 2 - 2A + £R + aa &*), (2.5) 

where R is the spatial Ricci scalar, a; = d r In N and A, A, £ and a play the role of the 
independent coupling constants we have mentioned. A is identified as the cosmological 
constant by analogy to GR. In (2.5) we have fixed the globally multiplicative constant of 
the action equal to one, since such a constant does not affect the vacuum held equations. 
For a similar reason we may fix the value of one more coupling constant since, under the 
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conditions we are going to impose on the field equations, one constant can be absorbed 
with rescalings of the other ones. 

Now we start with studying the solutions of the action (2.5) we are interested in. By 
the standard procedure of setting spherical coordinates and adapting the time coordinate 
to the timelike Killing vector, the most general static spherically symmetric spacetime 
metric can be written in time x spherical coordinates as 


ds 2 


~(N(r) 2 


N r N r )dt 2 + 2 N r drdt + 


dr 2 

f{r) 


+ r 2 dVt 2 2 ). 


( 2 . 6 ) 


N r = N r (r) is the component of the shift function that is left undetermined when only 
the conditions of staticity and spherical symmetry are imposed and not coordinate trans¬ 
formation has been yet invoked. As it is well known, in GR N r can be absorbed by 
a coordinate transformation that implies mixing time with space. However, in Horava 
theory such a transformation is not valid since it implies changing the foliation of the 
spacetime. Thus, the shift function remains as an active functional degree of freedom 
that distinguishes between inequivalent configurations among the set of static spherically 
symmetric configurations. For the sake of simplicity we shall set the shift function equal 
to zero as part of our ansatz, thus reducing the space of solutions. Note from (2.3) that 
staticity and vanishing of the shift function imply K t] = 0, and as a consequence of this 
A disappears from the field equations. Hence all the analysis is valid for the A = 1/3 and 
A 7 ^ 1/3 theories. The vanishing of all the kinetic terms leads to a further simplification 
in the set of coupling constant: one of the three remaining constants A, £ or a, can be 
absorbed by rescaling the other two, as we anticipated. From now on we set £ = 1. 
Therefore, we have that the vacuum field equations depend on two coupling constants, 
A and a , under the conditions we have set. 

From the action (2.5) we derive the field equations and then evaluate them on static 
configurations with vanishing shift function. This yields the equations 

R ij - ^ g ij R + A g ij - N~\V L V j N - g ij W 2 N) 

+aN~ 2 (y i N'V j N — ^g ij 'V k N'V k N) = 0, (2.7) 

R — 2A — a(2N~ 1 'V 2 N — N~ 2 X7iNV L N) = 0. (2.8) 


If a = 0 these equations coincide with the Einstein field equations with cosmological 
constant in the ADM formulation under the conditions of staticity and N t — 0. Equation 
(2.8) and the trace of (2.7) imply 

(2-a)N- 1 V 2 N + 2A = 0. (2.9) 


Thus, field equations admit no solution for a = 2 whenever A / 0. Assuming a ^ 2 , 
equations (2.8) and (2.9) are equivalent to 

V 2 N + ^N = 

R + aN~ 2 V k NV k N - (2 - 3a)y = 


6 


0, 

0, 


( 2 . 10 ) 

( 2 . 11 ) 



where 7 = 2A/(2 — a). After using these equations into the held equation (2.7), it 
becomes 

R ij - N^V^N + aN~ 2 YNV j N - (1 - a ) 7 g ij = 0 . (2.12) 

Thus, the system to be solved for static conhgurations with vanishing shift function is 
Eq. (2.12) together with (2.10) or (2.11). 

As we have discussed the general static spherically symmetric metric with a vanishing 
shift function is given by 


N = N(r), ds} 3) = j— + r 2 d(l 2 m . (2.13) 

Under this ansatz Eqs. (2.10) and (2.11) yields, respectively, 

,— N 

(<- 2 V7V)' +l-fj = 0, (2.14) 

rf + f - 1 - \r*S ((A + (1 - 3a/2) 7 r 2 = 0. (2.15) 

All off-diagonal components of the equation of motion (2.12) vanish. The rr and 66 
components become, respectively, 

0 , (2.16) 

0. (2.17) 


/' N" f'N' 

— -- 1 - -— cy. 

rf N 2fN 


N'y , . 1 

w) +(1 -“ h 7 = 


1 f , , rfN' 2 

2 r f + / - 1 + + ( 1 ~ a > rr = 


The system of equations (2.14 - 2.17) is redundant. Indeed, after doing some computa¬ 
tions, it turns out that the first-order equations (2.15) and (2.17) imply the second-order 
ones (2.14) and (2.16). Thus, the number of independent equations matches with the 
number of unknowns, N(r) and /(r). Since Eqs. (2.15) and (2.17) are first-order equa¬ 
tions there are two independent integration constants in the full set of solutions. One of 
these constants is evidently associated to scalings of N, which can always be absorbed 
by scaling the time. Therefore, only one integration constant has physical meaning. Al¬ 
though Eqs. (2.15) and (2.17) determine the closed system of equations, we shall strongly 
use the Eq. (2.14). 

The ansatz (2.13) includes the Schwarzschild-de Sitter and Schwarzschild-anti-de Sit¬ 
ter spacetimes. They are given by 


N 2 = f 



(2.18) 


where the Schwarzschild-de Sitter metric holds when A > 0 and Schwarzschild-anti-de 
Sitter for A < 0, and M is a free constant. If we set the GR value a = 0, which implies 
7 = A, these metrics solve exactly the Eqs. (2.15) and (2.17) for any A and M. 
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2.2 Deformed AdS asymptotia 

Before entering in the full solutions, we study the potential asymptotic divergences for 
the case when the cosmological constant is negative. This helps to contrast with the 
anti-de Sitter metric, which is (2.18) with A < 0 and M — 0. To achieve this we assume 
that if N and / diverge at r —y oo, then they do so with some dominant powers of r. 
That is, we assume that as r —» oo, 

N = r a , / = Cr b , (2.19) 

with a,b, C > 0 (again, Eqs. (2.15) and (2.17) give no information about multiplicative 
constants of N). This assumption is supported by the numerical solution, as we shall 
see. We then evaluate the field equations at the limit r —* oo by substituting (2.19) into 
them and neglecting any finite term. Equation (2.14) yields 



Ca{ 1 + a + b/2)r b + yr 2 = 0 . 

(2.20) 

Since none of these coefficients can be put equal to zero, we 

have that this equation 

necessarily implies 

b = 2 and Ca( 2 + a) = —7 . 

(2.21) 

Next, Eq. (2.17) yields 

C(1 + a + b/2)r b + 7(1 — a)r 2 = 0 . 

(2.22) 

With the same argument 

we have again that b = 2 and 



C{ 2 + a) = —7(1 — a). 

(2.23) 


We can solve a and C from (2.21) and (2.23), which yields 

a = ——— , C=--(- - ^ -— V (2.24) 

Finally, it can be checked that with the values for a,b,C given in (2.21) and (2.24) the 
Eq. (2.15) is solved asymptotically. Thus, we see that the solutions of the field equations 
that diverge asymptotically with some powers of r do not tend exactly to the anti-de 
Sitter space, but rather to a deformation of it (for small a). Specifically, / does diverge 
as in AdS, / ~ r 2 , but N exhibits a slight modification, N 2 ~ r 2(i-«) ^ which for 

small a can be seen as a deformation of r 2 . The field equations fix the coefficient of the 
dominant mode of /, which is the constant C given in (2.24) (actually, this coefficient is 
also different to the one of AdS). Thus, the dominant mode of / at infinity is completely 
fixed for all solutions 

2.3 The coordinate transformation 

Let us start our approach for solving the system of equations. We remark that in the 
whole procedure we shall assume that a is near to zero. There is a quadratic structure 




encoded in Eqs. (2.15) and (2.17). Indeed, a combination of these two equations yields 
the quadratic equation 




i, 


(2.25) 


where 

p 2 = 1 - Ar 2 , fj = y/l - a/2 . (2.26) 

Our assumption for a implies that j3 is close to one. This quadratic equation suggests 
the implementation of a procedure similar to the one we used in the A = 0 case [21]. 

To further proceed we choose a negative cosmological constant from now on. With 
A < 0 we have that p 2 > 0 everywhere, so the analysis of Eq. (2.25) simplifies greatly. 
Under this condition we may write the Eq. (2.25) as 


(VJ r -VJN' \ 2 _ ( 0ry/jN' \ 2 = 

V P pN ) V pN j 

The general solution to this equation can be written as 


PryfiN' 

pN 


shill x, 


V? fy/jN' 
P PN 


± cosh x i 


(2.27) 


(2.28) 


for any x G (—oo, +oo). Thus, we have that in principle there are two sets of solutions, 
each one corresponding to a choice of sign in (2.28). It turns out that the two sets of 
geometries are equivalent after coordinate transformations. We shall further comment 
on this later on. Thus, we may take (2.28) with the upper sign as the general solution of 
Eq. (2.27) without loss of generality. Solution (2.28) is then equivalent to 


P 

rN' 

~w 


cosh y — (3 1 sinh x , 


1 / cosh x 1 \ 1 

~P V sin li X PJ 


(2.29) 

(2.30) 


Now our aim is to use x as a new radial coordinate, obtaining, Erst of all, an equation 
for the coordinate transformation between r and x an d then equations for N and / in 
terms of X- For this goal we join Eq. (2.14) with the Eqs. (2.29) and (2.30) to form a 
system of equations for the functions r(x), IV(x) and f(x)- Notice that in this procedure 
we are replacing the original held equations (2.15) and (2.17) with the Eqs. (2.14) and 
(2.27), which form an equivalent system of two held equations. 

We hrst rewrite the differential equation (2.14) as an equation with x as the indepen¬ 
dent variable, obtaining 


d 

dx 



r 2 N dr 

+7 T7T 


0. 


(2.31) 
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This equation can be brought to the form 


/ 1 dN\ 
\N~dx ) 


r 2 p 



f 1 dN\ 
\N~ik ) 



and finally to the form 


1 dN 
N dr 



r 2 dr 
P dx 


= 0, 
(2.32) 



r dN 
N~ik ) P 



r dN 
N dr 

v7 

p 


(1 — 2A r 2 ) + yr 2 — 

d\ 

± r fvn t^dN 

d X [\ P ) \N dr 


rp 


(2.33) 


By substituting Eqs. (2.29) and (2.30) into this last equation and after a bit of algebra 
we obtain a differential equation only for r(y), 


dr f 1 coshy\ rp 2 
dx VP sinliy ) H 


where 


H = 1 - A 



cosh x 
/3 sinh x 


(2.35) 


Equation (2.34) determines the coordinate transformation from y to r. 

The space-time metric can be expressed in the y coordinate by using Eqs. (2.29) and 
(2.34), which yields 


ds 2 = -N{xfdt 2 + h(x)dx 2 + r(x) 2 dXl 2 2) , (2.36) 


where 


Kx) = 


r 2 p 2 


sinh 2 y H 2 


(2.37) 


Thus, the radial component h( y) is algebraically determined once r(y) has been obtained. 
Now we use the chain rule in Eq. (2.30) and substitute Eq. (2.34) into it, obtaining 


1 dN _ p 2 

nHx ~ ~W 


(2.38) 


This equation allows to find N{ y) once r(y) has been found from (2.34). Under this 
approach we have that the integration of Eqs. (2.34) and (2.38) leads to the two inde¬ 
pendent integration constants we discussed above. Finally, we have that the solutions of 
the field equations are determined by Eqs. (2.34), (2.37) and (2.38). 

We recall the reader that, to arrive at this point, in Eqs. (2.28) we chose the solution 
with plus sign. When the minus sign is taken the three equations determining the solu¬ 
tions are modified, but the reparameterization y' = —y brings them exactly to the form 
in (2.34), (2.38) and (2.37). Therefore, any choice of sign in Eqs. (2.28) leads to the same 
configurations modulo coordinate transformations. 
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We comment that at the point A = 0 Eq. (2.34) becomes 


dr 



cosh x 
sinh x 


(2.39) 


and its integral is 


r(x) = 


ke x IP 


(2.40) 


sinh x 

This is the coordinate transformation we found for the case without cosmological constant 

[ 21 ]- 


2.4 Analysis of the field equations 

Having found the differential equations that determines r(x), iV(x) and h(x), the rest of 
our procedure will consist of in doing analysis to extract special points of these equations 
and then performing numerical integration on them. For the integration the key equation 
is (2.34), since it is a single equation for r(x) and the solutions for N and h depend on 
it. In turn, for the numerical solution of Eq. (2.34) it is mandatory to analyze the several 
singularities this equation has as well as its critical points. From now on we assume a 
positive a, such that a ~ 0 + and j3 ~ 1 _ . In the Appendix we summarize the results for 
negative a. 

We start with the singularities. It is easy to see that at both asymptotic limits, 
X = ±oo, Eq. (2.34) implies that dr/d\ is positive, hence r becomes monotonically 
increasing at both ends of x- Thus, there is a divergence at x = +°° that yields r — oo, 
the spatial infinity. On the other hand, as x —>■ — oo the radius r necessarily tends to zero. 
We shall check this behavior at both asymptotic limits it in the numerical integration. 
We remark that at the divergence x = +oo the factor H diverges. 

Further divergences of dr/d\ arise for finite x (with no analogs in the A = 0 case). 
The key feature to localize these divergences is that they are just the zeroes of H. Indeed, 
another potential source of divergences is the factor 


cosh x 
sinh x 


(2.41) 


which diverges at x = 0, but it cancels out with the divergence of H at the same point. 
Note that this behavior is completely different to the singularity of the “Minkowskian” 
equation (2.39), which effectively diverges at x = 0. Therefore, characterizing the singu¬ 
larities of Eq. (2.34) for finite x is equivalent to study the zeroes H. 

We find that there are two kinds of zeroes of H depending on whether r diverges or 
not: 


1. Singularities of the first kind These are zeroes of H at which r diverges. Let 
us denote by Xs^ a zero of H of this kind. It turns out that the factor 



cosh x 
/3 sinh x 


(2.42) 
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arising in H must necessarily vanish at y., 1 ^ i n order to compensate the divergence 
of the r 2 factor that multiplies it, such that H remains finite. Therefore, a zero y^ 
of H is necessarily a zero of the factor (2.42). Since this factor does not depend on 
r its roots are fixed for all solutions, 

tanh y W = . (2.43) 

This has only one solution for any given j3, so there is only one y s 1 ' ) . With fixed 
we mean that this value is independent of any initial condition corresponding to 
any particular solution r(y); all the solutions r(y) that have this singularity reach 
it at the same point (2.43). Note that the converse for this behavior is not true 
in general: the factor (2.42) can be zero, but r may remain finite and in this case 
necessarily H = 1, so there is no divergence. That is, in the full set of solutions 
reaching the point (2.43), not necessarily all of them develop a divergence, some 
may pass through this point regularly. 

Since in this kind of divergence one has r = oo, it corresponds to the spatial infinity. 


2. Singularities of the second kind In the other kind of zeroes of H , denoted by 
X's \ r remains finite. In this case the factor (2.42) cannot vanish. Instead, one 
may obtain the value of r from Eq (2.35), 


- A(<f > Y = 


( cosh yi^ 


-l 




f3 sinh x. 


( 2 ) 


- 2 


(2.44) 


The value of x^ is not fixed in the set of solutions and the pair (xi 2 \r4 2 ^) corre¬ 
sponding to this kind of zeroes varies among each solution. Since relation (2.44) 
Exes the value of r in terms of the one of y, solutions passing through some y« ; 
may or may not exhibit the singularity depending on whether the image r(y.s ) 
coincides or not with rq In the last case H is nonvanishing. Thus, we have again 
that some solutions can posses this singularities whereas others do not. 

It comes out the question whether yi^ and y^ are just failures of the y coordinate, 
hence coordinate singularities, or essential singularities. In the case of y^ the asymptotic 
limit r —)■ oo is reached, hence there is a place for having at y —> y^ the deformed AdS 
asymptotia we have discussed. If this is the case yi^ would be a coordinate singularity 
and the divergence of the coordinate transformation (2.34) at yi^ would signal that we 
are reaching at one boundary of the space that is regular. In the case of y A s ’ we have 
that r acquires a finite value. If yi^ is also a coordinate singularity, such that the final 

/Q\ /Q\ /Q\ 

solution is regular at ys , then the solution written in terms of y will stop at rq yy^ ; ) 
but it could be continued if there are other available solutions ending regularly at the 
same value of rl; . We shall investigate these issues with the numerical evaluation of the 
metric components and their derivatives as well as the curvature. 

So far we have seen that the spatial infinity r = oo is reached in the y coordinate at 
X = Too and y = yi 1 . Are there other spatial infinities in the coordinate y?. We may 
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find analytically the asymptotic solution to the coordinate transformation (2.34). With 
it we may proof rigorously that the spatial infinity, r = oo, is reached in the coordinate 
X only at x = +oo and x = X^- This information is crucial to understand the behavior 
of the numerical solutions. To achieve this we first rewrite the Eq. (2.34) in the form 

dr = (3~ l - tanhx (2 

d\ 1 — Ar 2 p -2 (1 —/3 _1 tanh x) ’ 

and then consider the asymptotic limit r —y oo into it. p 2 can be substituted by p 2 = 
—Ar 2 , such that in the asymptotic limit the Eq. (2.45) becomes 


dr f3 1 — tanh x 
d\ 2 — /3 _1 tanhx 

This equation can be integrated directly. The result is 

k( 1 + tanh x) m 

r =- m -r , 

| tanh Xs — tanh x | n (1 — tanh x ) n 


(2.46) 


(2.47) 


where k is an integration constant and 


1 + (3 2f3 2 — 1 „ 1-/3 

— - 77 = - 77 = - 

4/3 + 2 ’ 4/3 2 — 1 ’ 4/3 - 2 ' 


(2.48) 


For j3 ~ 1 (a ~ 0 + ) we have m,n,n > 0. Therefore, Eq. (2.47) implies that r = oo 
is reached only at x = Xs^ and X = +oo. Moreover, one may estimate the error of 
the asymptotic solution (2.47) by comparing the equation it satisfies with the original 
equation. At r — y oo, the right hand sides of Eqs. (2.45) and (2.46) differ by terms of 
order 1/r 3 and lower. Thus, the error of the asymptotic solution (2.47) is of order 1/r 3 . 

Now we move to the critical points of the solutions of Eq. (2.34). For general (negative) 
A, equation (2.34) implies that r(x) has a finite critical point x a t the point where the 
factor (2.41) vanishes; that is, at 

tanh x — (3. (2.49) 

Since this condition is independent of r, we have that this critical point is also fixed. 
Surprisingly, this is the same location of the throat in the corresponding x coordinate for 
the A = 0 case [21], In the numerical solution we are going to see that x is effectively a 
critical point since r remains hnite (r/r is hnite) and H does not vanish at X- As for the 
X.s'^ and xi 2) singularities, there remains to elucidate whether x is a regular point of the 
solution, such that x just labels a regular point at which the r coordinate fails. This will 
be checked with the numerical evaluation of the metric components and their derivatives 
in the x coordinate. Taking a derivative to Eq. (2.34) and evaluating at x we get 


(l — 7(1 — a)r 2 ) 



sinh 2 x(r — Ar 3 ) . 


(2.50) 
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Recalling that we are considering a small a and a negative A, Eq. (2.50) implies that the 
second derivative of r(y) is positive at y. Therefore, y is a minimum of r(y). This shows 
that any solution possessing the point y l ias a throat at y. 6 

As a final analysis that will help to understand the numerical results, we study the 
ordering of these special points. First, by reversing the Eq. (2.44), we may write it as 

tanh x ( s 2) = (1-tanhyCb . (2.51) 

V 2Ar s / 

Since A < 0 and yi'" 1 > 0, this relation implies 0 < yi 2 -* < xi^ for all hnite r'p . We 

( 2 ) 

stress that this ordering holds regardless of the fact that the value of Xs is different for 
each solution whereas yi^ is hxed. Thus, yi^ plays the role of upper bound for all the 
possible values of y A s . Second, we recall that we are assuming a near to zero, which 
yields [5 close to one. Under this condition we have that y > yi 1 '* (this breaks down at 
/3 = l/\/2, far enough from /3 — 1). Therefore, we find the universal ordering 

0< xi 2) <xi 1} < X, (2.52) 

for a close to zero. Notice that no singularities are developed in the domain (y« ,+oo), 
except for the asymptotic singularity y —> +oo. 

2.5 Numerical solutions 

2.5.1 Global numerical characterization 

Now we are ready to perform the numerical integration. Equation (2.34) requires the 
specification of one initial condition r(y 0 ), where y 0 can take any value in (—oo,+oo). 
Since (2.34) with a given r(y 0 ) constitutes a closed initial-value problem, different curves 
r(y) cannot intersect themselves. We give several initial conditions by varying y 0 over 
the subsets specified in (2.52). In Fig. 1 we present the several kinds of solutions the 
numerical evaluation of Eq. (2.34) yields. 

There are three kinds of curves: 

1. Class 1: solutions with a minimum The first class of curves possesses the 
critical point y, whose location is the same for all this class of solution. They 
can be generated by giving initial data at the right of the first-kind singularity, 
Ao > As 1 '*- As we discussed above, there is no singularity in (yi^+oo), hence all 
these curves pass over y and get a minimum there. For y —>■ +oo they grow to 
infinity since they tend to the asymptotic singularity of Eq. (2.34). We see that at 
the left of y they pass through yi 1 without developing any divergence, but further 
at the left they meet a singularity of the second kind, such that the solution is valid 
up to this point and r remains hnite there. 

6 Note that this conclusion is based only on the unavoidable presence of a minimum of r( y) at the 
fixed value (2.49), with no assumption on the value of r itself, unlike the behavior of the singular points 
yi 1 "* and yi 2 ' 1 . Thus, any solution passing through y necessarily possesses a throat. 
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Figure 1: The several numerical solutions r(x) of the Eq. (2.34). The dynamical constants 
have been set as /3 = 0.9 and A = — 1, which yields the minimum approximately at x = 1-47 and 
the first-kind singularity at xi^ = 0.67. The solutions that develop the second-kind singularity 
Xs J are the blue and red curves. The Xs ; values, which the reader can locate at the left end 
of each of these curves, vary among each curve and all of them fall in the subdomain (0, xi^). 
There are three classes of continuous solutions, each one with a definite domain in x- The curves 
of class 1 (blue) start at their second kind singularity, pass through the critical point x, where 
they acquire a minimum, and then grow monotonically up to spatial infinity. The domain for 
these solutions is x 6 (x« ,+oo) and their images are composed of two sides. The image of 
X £ (x,+oo) is the side r E (r, oo) whereas the image of x £ (Xs ,X) is the side r E (f,r K s '), 
where r is the image of x- The class 2 (red) start at their second kind singularity and grow 
monotonically up to the first kind singularity where they reach the spatial infinity. They cover 
the domain x £ (xi 2 \xi^) and range r E (r.g 2 \ oc). The curves of class 3 (golden) are also 
monotonically increasing. They cover all the domain x £ (—oo^i 1 ^) with range r E (0, oo). 
Since 0 < xi^ < Xs^ < X-, the classes 1 and 2 never pass to negative Xi the curves of class 2 are 
always located in subdomains of (0, Xs ) an( 4 the classes 2 and 3 do not develop the minimum 
X- 


The numerical integration shown in Fig. 1 indicates that these solutions have two 
sides joined by a throat. In one side (from x — +°° to x = x) they cover the space 
from infinity down to the throat; in the other side they cover the space increasing 
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from the throat up to the finite value of r determined by the singularity of the 
second kind yj 2 . 

2. Class 2: intermediate curves The second class of solutions have both the sin¬ 
gularities of the first and second kind. They can be seen in the plot as the inter¬ 
mediate curves that start at some finite value of y located in (0, y^) and then 
increase monotonically. The points where they start are singularities of the second 
kind, so r remains finite there. These solutions are not valid for values of y lower 
than their corresponding Xs ; singularities. As y increases from the corresponding 
Xs ; the solutions grow monotonically up until they reach the singularity of the first 
kind y^\ which is universal and r = oo is reached. Thus, these solutions cover a 
subset of the space from the infinity down to the finite value rq ; determined by the 
singularity of the second kind they reach to the left. 

3. Class 3: Solutions with the origin The third class of solutions are the ones 
located at the left of the graph. They start at y —» — oo, where r = 0, and grow up 
until they reach the singularity of the first kind at yi 1 ^, where r = oo. They can 
be generated by given initial data at y 0 = 0. These solutions cover a single copy 
of all the possible space in r. They never develop a singularity of the second kind 
despite of the fact that they pass over values of y where solutions 1 and 2 reach 
these singularities. 

To obtain the metric components we must solve the Eq. (2.38) numerically. This 
requires the specification of a further (and last) initial condition, 7V(yo). In Fig. 2 we 
show the plots for the metric components N( y) and h( y) for the classes of solutions 1 and 
2 and in Fig. 3 for the class 3. Since the Eq. (2.38) for N( y) depends on r(y), To a given 
initial condition N(xo) there correspond several solutions (distinguished among them by 
r(y 0 )), hence the curves N( y) intersect themselves in general. The same happens with 
h( y) since its formula (2.37) is a noninjective function of r. For the sake of clearness we 
have plotted only a representative group of (almost) nonintersecting curves. 

The main features are: 

1. Class 1 As we have discussed the domain of these solutions is y e (y$ , +oo). The 
numerical integration shows that both functions N and h are completely finite and 
nonzero on this domain. In particular they are finite and nonzero at the minimum 
y, which in the plots is approximately at y = 1.47. This shows that the minimum 
y is a regular point where the coordinate r fails, since it repeats its values around 
y. Thus, there is a regular throat at y. Unlike r, The coordinate y is useful to 
parameterize the space around the throat. Since N and h are regular and nonzero 
in the inner domain (y^ 2 \ +oo), we conclude that there are no horizons or essential 
singularities in this domain (we recall the reader that the transformation from r to 
y determined by Eq. (2.34) was purely on the radial direction, with no dependence 
in time involved). Thus, any geodesic passes from one side of the throat to the 
other one without reaching horizons or essential singularities at the inner points. 
There is a minimal-area two-sphere of r = constant at the throat and the area of 
the spheres grows at the two sides of the throat. 
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Figure 2: Metric components In N and h for the classes of solutions 1 and 2. Panels (a) and (b) 

correspond to the class 1 and panels (c) and (d) to the class 2. In (a) and (b) each curve starts 

( 2 ) 

from the left at its corresponding x$ ; singularity, passes through the throat x and continue 

( 2 ) 

towards x +oo. N and h are completely regular in the whole domain x £ (Xs > +oo) where 
these solutions are defined, which includes the throat, h diverges as x goes to x.s , whereas N 
remains finite there. At x —> +oo N and h go to zero. In (c) and (d) each curve starts from the 
left at its corresponding xi^ singularity and continues towards the fixed x.s' "* singularity, which 
is located at the right end of all curves. N and h are regular in the full domain x £ (xi 2 \ Xs^)- 
h diverges at both singularities xi^ an d xi 2 ^ whereas N diverges only at xi^- The values of j3 
and A are the same of Fig. 1. 


For x higher than the throat the chart extends itself covering one side np to the 
spatial infinity r = oo (x = +oo). From the plots in Fig. 2 (a),(b) we see that N 2 
and h go simultaneously to zero at this asymptotic limit. We shall see that the 
Kretschmann scalar diverges at this limit, showing that this is an asymptotic es- 
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Figure 3: Metric components IniV (a) and h (b) for the solution of class 3 which contains the 
origin. The functions N and h are regular in the domain x £ (—o^xi^) which is in one-to-one 
correspondence with r £ (0,oo). N diverges and h goes to zero at x ~^ —■oo ( r —y 0); whereas 
N and h diverge as x — > xi^ ( r ~^ oo). The values of (5 and A are the same of Fig. 1. 


sential singularity. The behavior of N 2 and h at this asymptotic limit is completely 
equivalent to what happens in the infinite boundary of the inner side of the A = 0 
solution [15, 21]. 

In the other side of the throat, x £ (Xs ,X), the solutions stop at the second-kind 

singularity. From Eq. (2.37) we see that h(x) diverges at Xs ; since H = 0 there, 

which is confirmed in the plots. We did not End any divergence or vanishing of N 

at xi 2) numerically. Since N behaves as a scalar under pure spatial transformations, 

a feasible way to elucidate whether x'P is a coordinate singularity is to go back to 

the original coordinate r, under which / -1 is the corresponding radial component of 

the metric. We evaluate the function / numerically using (2.29). We find that / is 

finite and nonzero at Xs which we shall show explicitly in a procedure of joining 

solutions. The regularness of N and / at Xs shows that Xs is a coordinate 

singularity and r is a good coordinate to cover this point (notice that the validity 

( 2 } 

of the coordinates goes in opposite ways between \ and Xs ). 

2. Class 2 The domain of these solutions is x £ (Xs 2 \xs^); N an d h are finite and 
nonzero on this domain. Thus, there are no horizons or naked singularities in inner 
points. These solutions cover regularly an open subset of the space from the spatial 
infinity r = oo (x = xi^) down to a finite point r = r ^ (x = X^)- Both N(x) 
and h(x) diverge at the right of the plots, where the first-kind singularity x^"* is 
reached. This leads us to ask ourselves whether this asymptotic divergence follows 
the deformed AdS asymptotia we discussed previously. 

At the left end, where the solutions stop at Xs , h diverges again whereas N remains 


18 












finite and nonzero. Thus, we have that the behavior of this solution at the second- 

kind singularity is very similar to the previous class. We find that f(x^) is finite 

( 2 ) 

and nonzero, hence for these solutions x.s is also a coordinate singularity. 

3. Class 3 The domain of these solutions is y e (—oo, yi 1 "*); N and h are finite 
and nonzero over this domain. These solutions cover the full space r 6 (0, oo) in 
a regular fashion, there are no horizons or essential singularities in inner points, 
and there are no critical points. At the limit r —> 0 N diverges and h goes to 
zero. The divergence of N leads us to think that this is an essential singularity, 
which we are going to confirm with the curvature. At the limit r —)• oo we have 
that N 2 , h | r=oc = oo, thus we shall look for deformed AdS asymptotia also in this 
solution. 


2.5.2 The joining of two solutions 

The fact that solutions of classes 1 and 2 end regularly below and above of some finite r 
suggests that both solutions can be joined to form a single wormhole geometry composed 
of two sides joined by a throat (the throat class 1 has), with its two sides extending 
themselves to spatial infinity, as the wormhole solution of the A = 0 case [15, 21]. 

To clearly show the joining, we generate solutions r(y) of classes 1 and 2 that tend 
to the same Xs ; singularity by giving initial data in the form 

r (xi 2) + S) = rf ) ±e. (2.53) 


For the case —e we obtain curves with the throat y (class 1) and for +e we obtain the 
intermediate curves that have both kind of singularities (class 2). For 5 and e sufficiently 
small we get that both curves tend to join themselves at y^ , such that the image of 
their union, including the value A 2 \ constitutes an extended, continuous, and doubly 
copied range of r. In Fig. 4 we show such continuous joinings of the functions r(y;). We 
generate several curves by varying yi 2 in the range (0, y,, ). We see that in all cases the 
coordinate r is regular and continuous at the joining point. 

We may also see continuity and smoothness in the metric components. The derivative 
dlnN/dr is taken from Eq. (2.30) and df /dr follows by taking a derivative on Eq. (2.29) 
with respect to r and combining with Eq. (2.34). This yields 


dj_ 

dr 


H 


2 sinh y ( ; d 1 coshy — sinhyj-2A (coshy — [5 1 sinhy) r. 


(2.54) 


Thus, Eqs. (2.30) and (2.54) give dlnN/dr and df /dr as functions of y. In Fig. 5 we plot 
the numerical solutions for A r (y), /(yj and their derivatives using again (2.53) combined 
with N(Xs ; + 5) = N 0 ± e, where N 0 is a fixed value (this is equivalent to adjust the 
integration constants of both solutions). We may see that N , / and their derivatives are 
finite, nonzero and completely continuous at the joining point. Thus we conclude that 
solutions of the classes 1 and 2 form a single solution that is regular at the joining point 
ri 2 \ We remark again that the passage from y to r is a pure spatial transformation 
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Figure 4: Curves of classes 1 (blue) and 2 (red) that tend to the same second-kind singularity 

( 2 ) ( 2 ) 

Xs ; j plotted near the singularity. For each joined curve the Xs singularity corresponds to the 

joining point of its red and blue parts. 


involving only the radial direction. Thus, there is no horizon at the joining point. We 
have already discussed that the two separated solutions of classes 1 a 2 have no horizons 
or essential singularities at finite points. 

The theory has two static spherically symmetric solutions with vanishing shift func¬ 
tion: (i) a wormhole solution whose geometry consists of two sides, nonsymmetric between 
them, joined by a throat and with the two sides extending themselves from the throat to 
the spatial infinity. This solution is completely regular in its inner points, there are no 
horizons or essential singularities at finite locations, (ii) A solution that covers a single 
copy of the space from r — 0 to r — oo and that is regular in r e (0, oo). There remains 
to analyze the asymptotic limits on the two sides of the wormhole and the behavior of 
the last solution at r = 0 and r = oo, as well as the inspection of the curvature at all of 
these points. 

2.5.3 Numerical asymptotia and curvature 

We contrast with the Schwarzschild-anti-de Sitter space and check the asymptotia de¬ 
scribed in the preliminary analysis. The Schwarzschild-anti-de Sitter space arises in our 
approach by putting (3 = 1 (a = 0) in Eqs. (2.34), (2.38) and (2.29). All the analysis of 
the singularities of Eq. (2.34) still holds with (3 = 1. Thus, there is a singularity of the 
first kind at tanhyi 1 - = 1/2 and singularities of the second kind in (0,y^). Again, two 
solutions ending at the same Xs ; can be smoothly joined to form single solution. When 
(3=1 the critical point x moves to infinity, so there is no throat at all, as expected. In¬ 
stead, for y —>■ Too Eq. (2.34) implies dr/dx —» 0 _ . This implies that the function r(y) 
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Figure 5: Functions (a) In N, (b) /, (c) dlnN/dr and (d) df/dr for solutions of classes 1 
(blue) and 2 (red) near the second-kind singularity. All of these functions remain continuous as 
functions of r. 


has a lower asymptote at x +oo. This asymptote is the location of the horizon, hence 

this solution in the x coordinate covers the exterior region of the Schwarzschild-anti-de 

Sitter space. This fact is corroborated in Eq. (2.47), since h = 0 when f3 = 1, meaning 

that r = oo is reached only at y = yi*' 1 in this case. The class 3 solution, which extends 

itself to x —^ —oo, has the essential singularity (the Schwarzschild-anti-de Sitter space 

with negative mass). These results are completely analogous to the way we recovered the 

Schwarzschild solution in the A = 0 case [ 21 ]. 

In Fig. 6 (a) we show the functions r(y) for a wormhole solution and for the (positive 

mass) Schwarzschild-anti-de Sitter space. Each solution has been joined at its respective 
( 2 ) 

Xs ; singularity. In general these solutions are very different since the wormhole covers two 
copies of a subset of the space, whereas the exterior black hole covers only one. Looking 
at the divergences we have found so far, it is clear that the side of the joined wormhole 
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Figure 6: (a) Function r(x) with a = 0 (lower curve) and a ^ 0 (upper curve) in the range 

X £ (Xs ',+ 00 ). The a / 0 solution has a minimum, hence it has the throat, whereas the 
a = 0 solution has a lower asymptote; it is the exterior of the Schwarzschild-anti-de Sitter black 
hole. In (b) and (c) we compare N and / of one side of the joined solution with the deformed 
AdS asymptotia. The curve in (b) is In (N(r)/No(r)). It goes strongly to a constant value as r 
increases. The curve in (c) is /(r)//o(r). It goes strongly to one. 


solution that has some similarity with the exterior black hole is the one containing the 
intermediate solution (class 2). In Figs. 6 (b), (c) we compare the asymptotic behavior of 
the functions N(r) and f(r) on this side for a particular class 2 solution with respect to 
the asymptotia found from (2.19) to (2.24); that is, we compare with No = r and 
/o = Cr 2 , where C is given in (2.24). In the plot (c) we see that the ratio ///o goes very 
strongly to one. In (b) the logarithm of the ratio, In N/N 0} also goes very strongly to a 
constant, which is the expected behavior. The value of this constant does not need to be 
zero since it is the coefficient of the dominant mode of N, which can be regarded as an 
integration constant (this freedom is not present in /). Therefore, the plots in Figs. (6) 
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( 2 ) 

Figure 7: Kretschamnn scalar for the wormhole solution. The indicated points are A = In Xs 
(the joining point), B = In x^ and C = lny. We see that K. is hnite and continuous at the 
joining point. The asymptotic limit of the side of the solution that is asymptotically deformed 
AdS can be seen in the red curve, for which 1C tends to a finite value as r —> oo (x —> Xs^)- The 
asymptotic limit in the other side can be seen in the blue curve, for which /C diverges as r —> oo 
(x —>■ Too) . Thus, there is curvature singularity at this asymptotic limit. In this solution /C is 
regular over the whole solution, except at the essential asymptotic singularity at x —> Too. 


confirm the asymptotic behavior N 2 ~ r 2 ( 1_ “) 1 and / ~ Cr 2 in one of the sides of the 
wormhole solution. For small a, which is the case we have considered throughout all the 
analysis, we interpret this asymptotia as a deformed AdS space. In the other side, as we 
have already pointed out, N 2 ,h —> 0 asymptotically. For the class 3 solution we have 
also checked that there is deformed AdS asymptotia as r —» oo. 

The regularness of the curvature of the wormhole solution at the joining point as 
well as the potential presence of essential singularities in both solutions can be inspected 
in the Kretschmann scalar /C = R lil , a pR ,il ' a ^. In the x coordinate it takes the form 

1 / d 2 N ldhdN\ 2 8 f dN dh\ 2 

N 2 h 2 \ dx 2 hdx dx ) r 2 N 2 h 2 y ciy dx) 

2 / d 2 r 1 dh dr \ 2 4 A cZr \ 2 

r 2 h 2 \ dx 2 hdxdx) r 4 h 2 1 \dx) 

where r, N and h are understood as the functions of x that determine the metric com¬ 
ponents in (2.36). In Fig. 7 we show the plot for the joined wormhole solution (classes 1 
and 2). 

The first feature we remark is that the scalar /C is hnite and continuous at the joining 
point. This confirms that this point is just a failure of the x coordinate. Second, the 
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scalar /C is regular everywhere except at y +oo. At x +°° the spatial infinity on 
the inner side is reached. We observe that /C diverges at this limit, thus we conclude that 
there is an essential singularity at this asymptotic limit. At the asymptotic limit of the 
other side /C tends to a finite value, which is the behavior expected for the deformed AdS 
asymptotia. In the solution of class 3 we checked that /C diverges at r = 0 (y = — oo), 
thus there is a naked essential singularity at the origin in this solution, whereas at the 
r —» oo limit (y —* y« ) /C tends to a finite value which corresponds to the deformed 
AdS asymptotia, similarly to the exterior side of the previous solution. 


3 Conclusions 

We have focused the problem of finding the static spherically symmetric solutions of the 
Horava theory [1, 2] with a negative cosmological constant. We have taken the lowest 
order effective action of the complete nonprojectable theory and we have switched off the 
shift function in the spacetime metric (with only the symmetry of foliation-preserving 
diffeomorphisms available the spherical symmetry does not imply staticity and these are 
not enough to have a vanishing shift function). We have assumed a small a, which is the 
coupling constant of the (Sin N) 2 term. This term is the unique deviation from GR in the 
effective action [7], thus we are interested in a small a. We have carried out a handling 
of the held equations similar to the one of Ref. [21], where the same problem was focused 
for the case of vanishing cosmological constant. As it happened in that paper, here we 
have found a minimum for the physical radius r when it is regarded as a function of a new 
coordinate that runs at both sides of the minimum. This corresponds to a wormhole-like 
geometry of two sides joined by a throat, the minimum being the throat. The two sides 
are nonsymmetric between them and each one extends itself from the throat to the spatial 
infinity. The wormhole solution is everywhere regular and without horizons, except for 
an asymptotic essential singularity in one of the sides 

Curiously (and different to the A = 0 case), the new coordinate, y, fails at a certain 
finite point of one side. We have identified this as a coordinate singularity. Thus, nothing 
special occurs at that point but a failure of the y coordinate. Instead, the physical radius 
r is regular and valid there. We supported this conclusion by examining the metric 
components numerically and the square Riemann tensor (Kretschmann scalar). 

In the wormhole solution of the A = 0 case [15, 21] the asymptotia at the two sides is 
different. One side is asymptotically flat whereas the other one is asymptotically singular. 
The singular side can be regarded as a kind of inner space. In the case we studied here 
we found an analogous nonsymmetric asymptotic behavior, but with some differences. 
Asymptotically, one wish to compare the exterior side with AdS, which has N 2 , f ~ r 2 
asymptotically. However, we found analytically that the divergence of this side is of the 
form N 2 ~ r 2 (i~ a ) , / ~ r 2 . This was confirmed by the numerical analysis. Since a is 
considered small, we interpret this result as a deformed-AdS asymptotia, characterized 
by a growing of N 2 faster than AdS. The other side, which we regard as an internal space, 
exhibits an asymptotic essential singularity, as in the A = 0 case. 

Interestingly, the asymptotic behavior of the regular side is very similar to the Lifshitz 


24 


scaling solution, which has been used in the context of holographic duals of nonrelativistic 
configurations. In the context of the Horava theory, in Ref. [30] it was shown that the 
Lifshitz scaling spacetime is a vacuum solution of the complete nonprojectable theory with 
cosmological constant. To make a comparison with our solution, the “bulk” coordinate 
of the Lifshitz spacetime (r in the notation of [30]) can be equated to the radial direction 
we have used here. Although the spatial submanifolds of r = constant are different in 
both cases, the asymptotic tt and rr components of our solutions in the nonsingular side 
are exactly equal to the ones of the Lifshitz spacetime. In particular, the exponent of 
our —N 2 function, which motivates the name of deformed AdS, is equal to the exponent 
of the analogous component in the Lifshitz spacetime, which leads to the anisotropic 
scaling. Thus, asymptotically our solutions acquire a kind of Lifshitz scaling. 

Whole wormhole solutions with negative cosmological constant have been studied 
in GR and its extensions in many contexts. For example, they have been found in 
combination with scalar fields [31], their quantization has been carried out [32] and they 
has been used in the AdS/CFT correspondence [33]. We hope that the solutions of the 
Horava theory can find similar applications. 

The side with deformed AdS asymptotia can be used as exterior solution of stars, 
as was done in Ref. [15] for the analogous solution in the A = 0 case. This consists of 
taking the deformed-AdS side of the wormhole and truncating it at some finite radius 
greater than the location of the throat. Since the distribution of matter is located from 
the origin to the radius of joining, the vacuum solution is not valid inside the joining 
surface. Thus, such solutions have neither throat nor internal side at all. As exterior 
solution, it is interesting to have an asymptotic behavior slightly different to the one of 

AdS. 

Returning to the whole wormhole solution, an important feature it has is that it is 
a vacuum solution. In GR it is well known that exotic matter is needed to maintain 
a traversable wormhole. But in the solution we have at hand this is not the case, its 
existence does not require the coupling to matter. Geodesics pass from one side to the 
other one through the throat without seeing any matter, neither exotic nor ordinary. 
Thus, no exotic matter is needed to maintain this solution. The same feature holds 
for the wormhole solution of the A = 0 case. This is conceptually similar to some 
traversable wormholes that have been found in the braneworld scenario [34, 35], where 
the field equations allow the space to maintain its throat by itself, with no need of matter 
violating energy conditions (in the braneworld it is the effect of the bulk geometry that 
maintains the throat on the brane). 

We have also found solutions with naked singularities, like in the A = 0 case. In 
particular, only such solutions arise in the case of negative a. These solutions cover a 
single copy of the r e (0, oo) space and also describe deformed AdS asymptotia. 

About the name that can be given to our two-side solution, we comment that the 
notion of a wormhole geometry is a global one, it is not sufficient to find a throat, which is 
a local feature, to conclude that some space can be classified as a wormhole. In particular 
there should no be physical barriers avoiding to pass from one side two the other one in 
the two directions through the throat. An study on these issues for some solutions was 
carried out in Refs. [36, 37]. It was found that some solutions of the Brans-Dicke theory 
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and the f(R) theory proposed as wormholes are actually global wormholes only under 
conditions for which the theory propagates ghosts. In some instances configurations 
named as wormholes are considered with diverse asymptotia, see [38, 39]. The numerical 
analysis we have performed here constitutes a global study of the configurations. We 
have found no barriers that would avoid the two-way movement of physical particles 
between the internal points of the two sides of the solutions (excluding the asymptotic 
singularity at one side). We also comment that the complete nonprojectable Horava 
theory at A = 1/3 has no ghosts (it has no extra degrees at all) [6] and for A ^ 1/3 one 
may require A > 1 in the linearized theory to safe the extra degree from becoming a ghost 
[2], The solutions we have found here are valid for any A. 

Finally, the stability of our solutions under small perturbations is an important aspect 
that deserves a detailed analysis. Such an investigation is beyond the scope of this paper. 
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Appendix: Negative a 

When a is small and negative, a ~ 0 _ , we have that f3 ~ 1 + . In this case all the 
analysis remains identical up to the Eq. (2.40). The singularities of Eq. (2.34) are 
very similar. There are first-kind singularities at tanh y,, 1 ' 1 = 1/2 j3 and second-kind 
singularities in (0,yi^). The interpretation of both singularities is the same: they are 
coordinate singularities. Thus, the two solutions ending at the same y$ ' can be joined 
to form a single solution. Therefore, there are two solutions: one formed by the joining 
of two solutions in (yi 2 \ +oo) and the other one in (—oo, yi^). 

The sector y G (—oo,ys ) is very similar to the a > 0 case. These solutions have 
neither second-kind singularities nor critical points. They have an essential singularity 
at the origin and cover r G (0, oo). We confirmed this by evaluating numerical solutions. 

The sector (y« ,Too) is very different to the a > 0 case. The derivative dr(y)/dy 
at y —» Too becomes negative. Thus, r becomes monotonically decreasing at y —* Too, 
such that the spatial infinity is not reached at y —> Too, but instead the solution goes 
to the origin. This is confirmed by Eq. (2.47), since now n < 0. Moreover, there is no 
possibility of critical point here (Eq. (2.49) has no solution). This implies that this part of 
the solution is monotonically decreasing in the full domain (yi 2 \ Too). Hence, its union 
with the other part covers the full range r G (0,oo). It has a naked essential singularity 
at the origin. We confirmed this result with numerical integration. 

Therefore, in the a ~ 0“ case all the static spherically symmetric solutions with 
vanishing shift function have naked singularities at the origin and cover a single copy of 
the full space r G (0, oo). They acquire deformed AdS asymptotia as r —>• oo. 
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